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SUMMARY 


The  Electrical  Geophysics  Laboratory  of  The  University  of  Texas 
performed  a series  of  measurements  of  the  electromagnetic  emissions 
from  merchant  ships  passing  through  the  Corpus  Christi  Ship  Channel 
near  Port  Aransas,  Texas.  Analysis  of  the  recorded  signals  shows 
that  they  correspond  closely  to  the  fields  excited  by  a horizontal  electric 
dipole  source  moving  with  the  ship.  Also  recorded  were  the  levels  of 
the  natural  electromagnetic  background.  Using  these  levels  the 
maximum  ranges  at  which  a matched  filter  could  extract  the  ship  emission 
signals  are  derived.  For  conditions  typical  of  the  open  ocean  detection 
ranges  of  up  to  2 km  are  obtained.  Also  considered  are  the  effect  of 
electromagnetic  noises  generated  by  ocean  surface  waves  moving  the 
conductive  sea  water  in  the  presence  of  the  geomagnetic  field.  > 


lit 


I. 


INTRODUCTION 


During  the  months  of  September  and  October  of  1973,  the 

Electrical  Geophysics  Laboratory  of  The  University  of  Texas  operated 

sensitive  electric  and  magnetic  field  recording  equipment  at  a shore 

site  adjacent  to  the  Corpus  Christi  Ship  Channel  near  Port  Aransas, 

Texas.  The  purpose  of  the  experiment  was  to  record  the  electric  and 

magnetic  field  emissions  from  passing  merchant  ships  in  the  frequency 
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range  extending  from  approximately  10  Hz  to  5 Hz.  In  addition, 
measurements  were  made  of  the  background  noises  that  would  serve 
to  limit  the  minimum  levels  at  which  the  emissions  could  be  identified. 
From  these  measurements  maximum  ranges  at  which  the  ships  could  be 
detected  by  sensing  their  electromagnetic  emissions  have  been  com- 
puted . 

One  additional  consideration  was  included  in  the  range  computa- 
tions. This  involves  the  theoretical  computation  of  the  electromagnetic 
noisos  that  would  be  generated  by  ocean  waves  moving  conductive  sea 
water  in  the  presence  of  the  geomagnetic  field.  This  was  done  in 
anticipation  of  the  fact  that  an  actual  detection  station  may  be  situated 
at  sea.  At  a sea  site  the  wave  noises  would  add  to  the  natural  back- 
ground noises  to  further  limit  the  maximum  range  of  detection.  An 
appendix  presenting  the  theoretical  derivation  of  the  wave  noise  fields 
has  been  included  with  this  report. 
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II.  THE  MEASUREMENTS 
A.  General 

The  location  of  the  field  site  adjacent  to  the  Corpus  Christi 
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Ship  Channel  is  shown  in  Fioure  1 . Figure  2 shows  a plan  view  of  the 
field  site  layout  including  the  sensor  positions  and  the  ship  channel 
margins . 

During  the  period  September  18  through  October  7 , a total 
of  55  data  runs  were  recorded.  Of  this  total  37  involved  the  measure- 
ment of  the  electromagnetic  emissions  of  passing  ships,  13  involved 
recordings  of  the  background  noise  levels  when  no  ships  were  in  the 
vicinity  of  the  sensors , 4 were  calibration  runs  and  1 involved  the 
recording  of  parallel  electric  field  sensing  lines  to  test  the  noise 
levels  of  the  electrodes.  All  of  these  runs  are  summarized  in  Table  1. 
The  runs  made  during  ship  passages  are  designated  with  the  name  of 
the  ship.  The  other  runs  are  designated  as  to  type. 

Each  data  run  whether  made  for  the  purpose  of  sensing 
the  emissions  from  a passing  ship  or  to  determine  the  level  of  natural 
background  activity,  consisted  of  recording  the  data  in  two  forms. 

One  was  a strip  chart  record  of  the  analog  signals.  The  other  was  a 
digital  tap*,  jcording  of  the  data  sampled  and  digitized  at  0.1  second 
intervals.  The  digital  tapes  contain  ten  minutes  of  data  giving  a 
total  of  6,000  points  per  run  for  each  component  signal  recorded. 
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During  ship  runs  the  ships  were  tracked  by  visual  sightings 
made  with  a surveyor's  transit.  The  location  of  the  transit  is  indicated 
in  Figure  2.  The  time  of  visual  passage  of  the  ship's  bow  or  stern, 
depending  upon  whether  the  slant  range  was  closing  or  opening 
respectively,  was  recorded  at  angles  of  30°  and  60°  either  side  of  a 
0°  or  perpendicular  sighting  on  the  channel.  The  times  of  both  bow 
and  stern  passages  were  recorded  at  the  0°  position. 

B . Equipment 

The  equipment  used  to  record  both  the  electromagnet’’ _ 
ship  emissions  and  the  background  noises  consisted  of  apparatus  to 
sense  both  the  electric  and  magnetic  fields.  This  apparatus  was 
designed  and  built  by  the  Electrical  Geophysics  Laboratory  and  has 
been  used  to  record  low  level  fluctuations  in  the  geomagnetic  field 
(geomagnetic  micropulsations)  and  the  associated  electric  fields 
resulting  from  the  induced  earth  currents  (telluric  currents)  in  a variety 
of  applications . 

Each  of  two  orthogonal  electric  field  components  tangent 
to  the  earth's  surface  is  sensed  indirectly  by  measuring  the  potential 
difference  between  earthed  electrodes  implanted  at  each  end  of  a 
baseline.  The  potential  difference,  V,  measured  between  the  two 
electrodes  is  related  to  the  electric  field  by  the  expression 


where  the  path  P follows  the  wires  connecting  the  two  electrodes  to 
the  potential  sensing  apparatus.  The  potential  sensing  equipment 

f 

consists  of  differential  input  amplifiers  in  cascade  with  a network  of 
post  amplifiers  and  filters , the  purpose  of  the  latter  being  to  amplify 
and  condition  the  measured  potential  difference  for  digitization  and 
storage . 

Each  of  three  orthogonal  components  of  the  magnetic  field 
are  sensed  with  one  of  three  induction  magnetometers  that  are  con- 
structed by  winding  a coil  of  copper  wire  around  a mumetal  bar.  The 
high  magnetic  permeability  of  the  mumetal  serves  to  concentrate  the 
flux  in  the  windings  and  increase  the  sensitivity  of  the  induction 
sensor.  The  aspect  ratio  of  the  bar  has  been  designed  so  that  its 
demagnetization  factor  desensitizes  the  sensor  to  variations  in  the 
permeability  of  the  mumetal  caused  by  various  influencing  environmental 
factors.  An  absolute  calibration  accuracy  of  2%  has  consistently  been 
observed  for  these  induction  magnetometers. 

The  voltages  induced  in  the  magnetic  field  sensors  are 
amplified  in  solid  state  chopper  type  preamplifiers.  These  amplifiers 
have  been  specially  designed  to  provide  very  low  noise  amplification 
of  the  induced  signals  at  the  impedance  levels  presented  by  the  sensors. 
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Post  amplification  and  filtering  is  similar  to  that  used  in  the  electric 
field  system. 

Ill . DATA  PRESENTATION  AND  ANALYSIS 
A.  General 

Of  the  total  of  37  runs  made  during  ship  passages,  the  10 
exhibiting  the  highest  quality  data  were  selected  for  analysis.  The 
quality  of  the  data  were  measured  by  the  absence  of  extraneous 
cultural  noises  or  irregularities  caused  by  equipment  malfunctions. 

No  judgement  was  made  on  the  basis  of  the  size  or  character  of  the 
emission  signals.  As  may  be  seen  in  the  data  presentation  to  follow 
some  of  the  ships  produced  large  emission  signals  while  for  others 
almost  none  were  detected.  One  of  the  most  persistent  sources  of 
extraneous  noise  was  the  presence  of  moving  automobiles  in  the 
vicinity  of  the  sensors  during  a ship  passage.  The  field  site  was 
located  in  an  area  of  free  public  access  and  it  was  considered  im- 
possible to  stop  all  vehicular  traffic  during  a data  run. 

Of  the  13  data  runs  made  on  the  background  noise  l*^vels 
10  were  considered  of  high  enough  quality  to  receive  further  analysis. 
The  results  of  only  4 of  the  10  background  runs  analyzed  are  presented 
here.  It  was  felt  that  this  number  was  adequate  to  characterize  the 
different  type  of  noise  fields  observed. 


B.  Physical  Ship  Parameters 

The  transit  sightings  made  during  ship  passages  were 
analyzed  to  yield  the  ship  parameters  length,  velocity  and  accelera- 
tion. For  the  analysis,  it  was  assumed  that  each  ship  followed  a 
course  lying  along  the  center  line  of  the  channel.  The  position,  x,  , 

D 

of  the  bow  of  the  ship  relative  to  its  position,  x , at  the  time  of  the 

o 

first  60°  sighting  with  closing  range  is  assumed  to  be  given  by 

a 2 

X = X + vt  + — t 
DO  I 

The  position,  x , of  the  stern  with  respect  to  the  position  x is 
s o 

expressed  by 

a 2 

X =X  - L + Vt+-r-t 
SO  2 

The  parameters  x^,  v (velocity  at  t=0),  a (acceleration)  and  L 
(ship's  length)  were  then  determined  by  solving  for  those  values  that 
produce  a minimum  mean  square  best  fit  to  the  observed  position  data. 
The  parameters  thus  determined  for  each  of  the  10  ships  for  which 
data  were  analyzed  are  presented  in  Table  2.  The  standard  deviation 
of  the  estimated  length  parameter,  L,  is  also  indicated  to  provide  a 
measure  of  the  scatter  in  the  data.  Also  included  in  the  tabulations 
of  Table  2 is  the  direction  of  the  ships  course. 
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C.  Spectral  Analysis 

The  first  step  in  the  spectral  analysis  of  the  recorded 
data  consisted  of  Fourier  transforming  the  sampled  time  series  for 
each  electromagnetic  field  component.  A Fast  Fourier  Transform 
routine  was  used  and  this  required  that  each  series  contain  a number 
of  sample  points  equal  to  an  even  power  of  2.  Since  each  data  run 
produced  series  that  contained  6,000  points  a decision  had  to  be 
made  whether  to  expand  each  series  to  8,192  points  by  the  addition 
of  zeros  or  to  truncate  the  6,000  available  points  to  4,096.  The 
second  alternative  was  adopted  because  it  was  reasoned  that  the 
dominant  spectral  components  would  be  generated  when  the  ship  was 
closest  to  the  sensors  and  this  period  of  time  was  relatively  short 
compared  to  the  total  length  of  time  required  for  the  entire  10  minute 
data  run. 

The  power  density  spectra  computed  from  the  Fourier 
transforms  of  the  component  fields  recorded  during  the  10  ship 
passages  and  3 background  runs  selected  for  presentation  are  shown 
in  Figures  16  through  28.  These  spectra  have  been  smoothed  to 
reduce  the  number  of  points  necessary  for  visual  presentation.  The 
smoothing  filter  has  an  approximately  constant  percentage  bandwidth 
of  26  percent. 
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D.  Strip  Chart  Recordings 

Reproductions  of  the  original  strip  charts  for  the  1 0 

selected  ship  passages  and  3 background  noise  runs  are  shown  in 

Figures  3 through  16.  Each  chart  includes  a coded  time  signal  for 

temporal  reference.  In  the  case  of  a shin  passage,  the  slant  range 

from  the  midpoint  of  the  sensor  to  the  ship’s  stern  is  indicated  for 

each  component  signal.  The  sensitivity  of  each  trace  is  also 

indicated  on  the  chart.  The  units  of  measure  are  those  commonly 

used  in  geophysics.  The  electric  field  intensity,  E,  is  measured 

in  millivolts  per  kilometer  (mv/km)  and  the  magnetic  flux  density,  B, 

-9 

is  measured  in  gamma  (Y)  where  ly  = 10  Webers  per  square 


meter. 


As  was  pointed  out  in  Section  IIB,  the  electric  field  was 


sensed  indirectly  by  measuring  the  potential  difference  between  two 
separated  electrodes.  It  is  this  potential  difference  that  is  actually 
recorded  on  the  strip  charts.  The  indicated  E field  units  have  been 
obtained  simply  by  dividing  the  measured  potential  difference  by  the 
straight  line  distance  between  the  sensing  electrodes.  This  is  an 
accurate  representation  of  the  electric  field  only  if  that  field  is  uniform 
over  the  span  of  the  array  and  if  there  is  negligible  extraneous  voltage 
Induced  in  the  electrode  wires.  This  extraneously  induced  voltage 
is  caused  by  the  time  rate  of  change  of  magnetic  flux  linking  the  loop 
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formed  by  the  straight  line  path  between  the  electrodes  and  the  actual 
path  of  the  wires  connected  to  the  electrodes.  For  the  slow  rate  of 
changes  of  the  magnetic  field  encountered  at  the  low  frequencies 
involved  in  this  study,  this  induced  voltage  is  negligible  for  the  wire 
layouts  actually  used. 

It  should  be  pointed  out  that  the  signals  recorded  on  the 
strip  charts  have  been  filtered  in  an  attempt  to  condition  (prewhiten) 
on  the  average  the  signals  for  digitization.  The  sensitivities  indicated 
on  the  charts  pertain  only  to  the  midband  regions  of  these  filters.  In 
some  Instances,  the  character  of  the  signals  outside  this  passband  has 
been  rather  severely  modified. 

E.  Analysis  of  Electromagnetic  Emissions 

Reference  to  the  strip  chart  reproductions  in  Figures  3,  4 
and  10  shows  that  the  emissions  from  the  three  ships.  Stainless  Trader, 
Benja  River  and  Mariotte  are  particularly  strong  and  that  most  of  the  com- 
ponent signals  can  be  distinguished  from  the  noise  to  slant  ranges  in 
excess  of  3,000  ft.  Because  these  signals  are  so  clear  it  was  decided 
to  compare  their  envelope  functions  with  those  that  would  be  produced  by 
a horizontal  electric  dipole  source  following  the  same  course  as  the  ship. 
Such  a source  is  chosen  for  comparison  "^ince  it  is  assumed  that  the  re- 
corded emissions  are  excited  by  the  ship's  propeller  modulating  the 
corrosion  currents  flowing  in  dipolar  fashion  between  different  parts  of 
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Dipole  Source 
ot  (0,0,0) 


Figure  III-l 


the  ship.  In  computing  the  signals  generated  by  the  dipole  source  it  is 
assumed  that  the  dipole  is  situated  at  the  water's  surface  at  the  position 
of  the  ship's  stern  and  that  it  is  oriented  parallel  to  the  keel.  These 
assumptions  represent  approximations  but  the  errors  that  they  cause 
should  rot  seriously  affect  the  principal  features  of  comparison  of  the 
computed  and  measured  signals. 

As  shown  in  Figure  III-l , the  dipole  of  moment  M follows 
a straight  line  path.  The  earth  is  assumed  to  be  uniformly  conducting 
with  conductivity,  a,  and  magnetic  permeability,  equal  to  that  of 
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free  space.  The  dipcie  provides  a quasi  sinusoidal  time  signal  with 
radian  angular  funuamental  fi'equency  At  the  point  of  closest  approach, 
it  passes  a distance  d from  an  observer  point  where  the  fields  are  computed, 
A rectangular  Cartesian  coordinate  system  with  axes  x,  y and  z is 
positioned  at  the  dipole  source  so  that  the  x axis  points  in  the  direction 
of  motion.  The  positive  z direction  is  up.  The  coordinate  system  is 
assumed  to  move  with  the  source  with  the  dipole  moment  pointing  in  the 
positive  x direction. 

Simulating  the  measurements, the  field  components  are  com- 
puted in  axes  parallel  ([j)  and  perpendicular  ( x)  to  the  dipole  path. 

As  may  be  seen  in  Figure  III-l,  the  )j  and  J-  axes  are  parallel  to  the 
X and  y axes,  respectively. 

Formulas  for  the  field  components  at  the  observer  point  are 
derived  from  expressions  presented  by  Bannister,  [Bannister,  1966]. 


The  derived  equations  are 


2 TTor 


(3  cos^  0-2)  + (1  + yr)  e J 


E = 


sin  0 cos  0 


ZTfor' 


|i  M 
o 


[«jK,  - 


(I  K, 
o 1 


- I.K  ) 1 sin  0 cos  0 (3) 

1 o J 


li  M 


S [(4  sln^  0 - 1)  I,Kj  - yf  (I^Kj  - IjK^)  sln^  ^] 


where 


[3  - (3  + 3 yr  + y^r^)  e 1 sin  0 (5) 

2 TT  Y r ‘- 


and  I , I, , K and  K,  are  modified  Bessel  Functions  with  argument 
o 1 o 1 

yr/2.  Equations  (1)  through  (5)  actually  represent  the  envelope 
functions  of  the  computed  signals  and  must  be  multiplied  by  a sinusoidal 
time  factor  to  simulate  the  recorded  quasi-sinusoidal  time  traces. 

The  computed  dipole  fields  are  compared  to  the  measured 
emission  signals  by  comparing  the  magnitudes  of  the  normalized  envelope 
functions . The  envelope  functions  for  the  data  are  scaled  from  the  strip 
chart  records  and  are  normalized  by  dividing  by  the  peak  value  during 
passage.  The  envelopes  of  the  computed  signals  are  also  normalized 
by  the  peak  value  attained  as  the  dipole  moves  past  the  observer  point. 

The  magnitudes  of  the  normalized  envelope  functions 
measured  from  the  three  ships  Stainless  Trader,  Benja  River  and  Mariotte, 
are  compared  with  the  same  functions  computed  for  the  moving  dipole 
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source  in  Figures  29,  30  and  31,  respectively.  The  experimental  data 


are  plotted  as  individual  points. 


3 

s 


The  theoretical  curves  are  displayed 
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with  a solid  line. 

The  fundamental  frequency  of  the  recorded  quasi-sinusoidal 
signals  was  som^-./hat  different  for  each  of  the  three  ships  and  also 
varied  to  some  extent  with  time  as  each  ship  passed  by.  The  interval 


over  which  the  frequency  varied  extended  from  1.2  Hz  to  2 Hz. 


For 


1 


the  computed  signals  a value  of  1.5  Hz  was  used.  A value  of  one  j 

mho  per  meter  was  used  for  the  conductivity,  a.  This  is  somewhat  | 

smaller  than  the  conductivity  of  the  sea  water  but  is  probably  near  the  j 

value  of  the  very  porous  sediments  both  underlying  and  adjacent  to  the  t 

ship  channel . I 

The  agreement  between  the  measured  and  computed  data  | 

is  seen  to  be  excellent.  The  small  differences  could  easily  be  caused  | 

by  the  ships  following  a course  slightly  to  one  side  or  the  other  of  the  | 

ship  channel  or  by  the  small  differences  in  the  fundamental  frequency  | 

of  the  emission  signals.  In  the  opinion  of  the  authors,  this  close  | 

agreement  between  theory  and  experiment  provides  clear  evidence  that  | 

the  signals  behave  like  those  of  a classical  horizontal  electric  dipole.  I 
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IV.  OCEAN  V^'AVE  NOISE  ANALYSIS 


A.  Electromagnetic  Noise  due  to  Ocean  'Waves 


Equations  for  the  electric  and  magnetic  fields  generated  by 
ocean  waves  in  the  conductive  sea  water  have  been  derived  using  simple 
Airy  surface  wave  theory.  The  derivation  is  presented  in  Appendix  A.  The 
Airy  theory  deals  with  a wave  for  which  the  water  surface  varies  sinusoidally 
with  both  time  and  space.  If  such  a wave  arbitrarily  propagates  in  the 
+ X horizontal  direction  and  the  water  motions  within  the  wave  are  indepen- 
dent of  the  horizontal  y direction,  then  the  water's  surface,  r\ , is  expressed 


j (ujt  - kx) 


Ti(x,t)  = ae 


where  g.  is  the  peak  amplitude  of  the  surface  undulation;  The  symbol  t 
represents  time  and 


is  the  propagation  constant  associated  with  the  wave  where  the  wavelength 


is  X.  For  this  study  the  water  is  assumed  to  be  deep  enough  to  allow  the 
"deep  water  approximation"  to  be  used  to  derive  the  dispersion  relation 


where  g is  the  acceleration  due  to  gravity. 
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The  subsequent  derivation  of  the  electromagnetic  fields  induced 
by  the  wave  motions  perpendicular  to  the  geomagnetic  field,  "5^  , separates  the 
induced  fields  into  two  parts.  The  sum  of  these  two  parts  constitutes  the 
whole  field  generated  by  the  ocean  wave.  One  of  the  parts  has  a magnetic 
field  that  is  transverse  to  the  direction  of  propagation  of  the  wave  and  to  the 
water  motions  within  the  wave.  For  the  wave  of  equation  (7)  this  means  that 
the  magnetic  field  has  only  a y component.  It  follows  that  the  fields 
associated  with  this  part  are  designated  Transverse  Magnetic  (TM) . The 
other  part  of  the  fields  produced  by  the  wave  of  equation  (7)  has  cnly  a y 
component  of  the  electric  field.  The  fields  associated  with  this  part  are 
called  Transverse  Electric  (TE) . The  field  equations  derived  in  Appendix  A 
for  these  fields  are: 

TM  fields: 


E = 

X 

j n 

juja  B e 
- oy 

(9) 

E = 
z 

n 

-u)a  B e 
- oy 

(10) 

B = 

V 

0 

(11) 

TE  fields: 
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I 


B = -ak(B  +JB  ) -/^,v 

X — ox  02  (k+jO 


(13) 


B = -jak  (B  + jB  ) [e^^  - ..  r 
2 — ox  02  (k  + jC) 


(14) 


where 


jC 


= y?7 


jl«UO 


(15) 


The  constant  jC  is  a complex  quantity  that  lies  in  the  first  quadrant.  The 
vertical  coordinate  2 is  positive  upwards.  It  should  be  pointed  out  that  the 
transverse  magnetic  flux  density  of  equation  (11)  is  2ero  as  a consequence 
of  neglecting  displacement  currents  in  favor  of  the  conduction  currents  in 
the  sea  water.  This  assumption  is  well  justified  at  the  low  frequencies 
involved  in  this  study . 

Equations  (9)  and  (10)  show  that  the  TM  electric  field  is 
circularly  polari2ed  and  is  always  parallel  to  a plane  of  propagation  of  the 
water  wave.  In  this  case  a plane  of  propagation  is  any  vertical  plane  con- 
taining the  direction  of  propagation  of  the  wave.  The  subscripts  n and  ^ are 
used  to  denote  field  components  that  lie  paralle  or  perpendicular  respectively 
to  a plane  of  propagation.  The  magnitude  of  the  rotating  TM  electric  field 
vector  is  given  by 

E||  = uua  B e^^  (16) 

The  TE  electric  field  described  by  equation  (7)  is  linearly  polari2ed  perpendi- 


cular to  the  plane  of  propagation.  Its  magnitude  is  expressed  by 


El  = uua  B 


ol 


I kz  2k  jCz 

I"  -CTTc®' 


(17) 


The  TE  magnetic  flux  density  like  the  TM  electric  field  is  circularly  polarized 
in  the  plane  of  propagation  and  its  magnitude  is  expressed  by 

B||  = akB  (18) 

" — on  ' k + ]C  ' 


Equations  (16),  (17)  and  (18)  form  a complete  description  of  the  magnitude  of 
the  electromagnetic  fields  generated  by  an  ocean  wave  having  the  sinusoidal 
behavior  considered  here.  Before  any  attempt  is  made  to  evaluate  the  wave 
induced  fields  for  typical  sea  states  we  must  first  inject  the  expressions  for 
the  full  spectrum  of  ocean  waves  that  are  present  in  the  real  ocean. 

B.  Ocean  Wave  Spectra 

A commonly  used  [ Conference,  1963  ] expression  for  the 


power  density  spectrum,  E(f),  of  ocean  wave  heights  is 

-n 


. -m  --[(f)  -1] 

E(f)  = E(y  (f ) e " o 
o 


(19) 


The  parameters  m and  n control  the  shape  of  the  density  function.  The 

variable  f is  frequency  in  Hertz  and  f^  is  the  frequency  of  the  mode  (maximum) 

of  the  density  function.  E(f^)  is  the  spectral  density  at  the  mode  frequency. 

It  should  be  pointed  out  here  that  E(f)  refers  to  the  spectral  density  of  peak 

2 

amplitudes  of  the  ocean  waves.  The  mean  square  or  variance,  a (f) , of  the 
ocean  waves  is  only  one  half  the  value  of  E(f): 


c^f)  = ^ 


(20) 
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Various  investigators  have  used  different  values  of  m and 
n to  match  individual  sots  of  experimental  data.  A sample  of  the  values 
reported  in  the  literature  [Conference,  1963]  is  shown  in  Table  3. 

Table  3 


Investigator 

m 

n 

Neumann 

6 

2 

Burling 

5.5 

7-9 

Roll-Fisher 

5 

2 

Bret  Schneider 

5-9 

4 

Although  the  choice  is  rather  arbitrary  the  values  used  by  Neumann 

will  be  used  here.  These  values  represent  the  spectrum  for  a fully 

arisen  wind  driven  sea . Substituting  m = 6 and  n ==  2 into  equation  (19) 

, -2 


gives 


, -b  " 

E(f)  = E(y(j— ) e 

o 


(21) 


The  mode  frequency  for  t;  e Neumann  spectrum  is  related  to  the  wind  speed,  v, 
through  the  expression  [Defant,  1961] 


o 2ttv  y 3 


(22) 


where  as  indicated  previously  g is  the  acceleration  due  to  gravity.  The 
spectral  density  at  the  mode  frequency  is  given  by 

-6 


E(f  ) = Cf 
o o 


(23) 
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where 

C = 1.37  X 10~^  (mVnz)  (24) 

gives  good  agreement  with  experimental  observations. 

We  may  derive  the  total  wave  amplitude,  E,  by  integrating 
the  density  function  to  obtain 

03 

E = J E(f)  df.  (25) 

o 

Substituting  equation  (21)  into  equation  (2  5)  and  integrating  gives 

E = f^E^(yK  (26) 

where 

“ 2I-  y 7 = -856  (27) 

The  most  important  characteristics  of  the  waves  can  be  found  in  terms  of  E 
as  [ Defant,  1961  ] 


most  frequent  value 

= 1.414 /E 

average  value 

H = 1.772  /E 

significant  value 

H = 2.832  JT 

^1/3 

average  of  highest  10% 

H = 3.600  /e" 

^1/10 

C.  Wave  Induced  Electromagnetic  Noise  Spectra 


Expressions  for  the  power  density  spectra  of  the  wave  induced 
noises  can  now  be  obtained  by  combining  the  wave  spectra  of  Sec  V-B  with 
the  electromagnetic  noise  expressions  in  equations  (11), (12)  and  (13).  In  so 


doing  we  note  that  the  induced  fields  have  been  derived  in  terms  of  a sinu- 
soidal water  wave  with  peak  amplitude  a.  As  derived  each  electro- 
magnetic expression  yields  the  peak  value  of  the  indicated  field  component 
where  that  component  field  varies  sinusoidally  in  response  to  the  water  wave. 
We  can  compute  the  power  spectral  density  of  the  wave  induced  field  compo- 
nents by  first  squaring  each  expression  and  then  replacing  the  square  of  the 
2 

wave  amplitude,  a by  the  appropriate  spectral  density  term.  For  our  purposes 

in  this  study  we  choose  to  compute  the  spectral  density  of  the  mean  square  or 

variance  of  the  v/ave  induced  fields  and  this  is  obtained  by  letting 

2 2,- 
a = a (f) 

Substituting  from  equation  (15)  for  a(f)  obtains 

= f (28) 

D . Wave  Induced  Electromagnetic  Noise  Computations 

The  expressions  for  the  power  density  of  the  wave  amplitude 

and  the  wave  induced  noise  spectra  were  programmed  for  a programmable 

calculator  and  the  spcctRim  values  computed  as  a function  of  frequency  for 

different  values  of  wind  speed.  Plots  of  the  computed  spectra  E(f),  P„  , 

^11 

P_  and  P_  are  shown  in  Fig.  32  through  Fig.  38  respectively.  Here 
Ej.  B|| 

the  symbol  P with  a subscript  is  used  to  designate  the  power  density 
spectrum  of  the  field  component  designated  by  the  subscript.  The  geophysical 
units  of  millivolts  per  kilometer  and  gamma  have  been  used  for  the  electric 
and  magnetic  fields  respectively. 
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.1^ 


The  terms  reprepjntinq  the  geomagnetic  field  B and  B 

oil  Ox 


are  each  a function  of  position  on  the  earth's  surface  as  well  as  the  direction 
of  propagation  of  the  water  waves.  Rather  than  compute  the  noise  fields  for 


a variety  of  values  for  B and  B it  was  decided  to  use  those  values  that 

oil  Ox 


would  produce  the  largest  noises  and  let  this  constitute  a "worst  case"  value. 


The  largest  value  for  B^^  would  exist  at  the  geomagnetic  poles  where 


B = 6 X 10  gamma 
oil 


(29) 


is  typical.  The  largest  value  for  B^^  would  occur  at  the  magnetic  equator 


for  waves  propagating  magnetically  east-west.  Typical  of  this  field  is  the 
value 


B^^  = 3 X 10  gamma 


(30) 


The  values  equations  (2^  and  (30)  were  used  for  the  computations. 

A summary  of  the  equations  used  to  compute  the  spectra  of 


Figs.  32  through  38  in  the  units  used  for  the  plots  is  given  below. 


(m  /km)' 


« in-6  2 E(f).„  2 2kz  ^ 

= 10  uj  Be  - — 

Ell  2 Ox  Hz 


(31) 


(m  /km)‘ 


Hz 


(32) 


B|| 

E(f) 


k+jC 


^2  Effl  ^2 

Hz 


-6 


f ,-2 


= E(f)(X, 


o' ' f 


Hz 


(33) 

(34) 


where: 
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V 


DETECTION  RANGE  COMPUTATIONS 


The  emission  signals  recorded  from  passing  merchant  ships  have 
been  shown  to  agree  closely  with  the  fields  of  a moving  horizontal  electric 
dipole.  This  means  that  the  dipole  fields  can  be  used  to  compute  the 
maximum  ranges  at  which  the  ship  emissions  could  be  detected.  For  the 
computations  it  is  assumed  that  the  ranges  between  the  ships  and  the  sensors 
would  be  somewhat  greater  than  that  which  existed  at  Port  Aransas.  For 
larger  values  of  y,  Equations  (1)  through  (6)  take  the  somewhat  simpler 
forms,  [Bannister,  1966]. 


M cos  w 


E = 


2 TT  ad 


^[(u-2)(u-n) 


(35) 


M cos  u)  t 


E = 

2 TT  ad" 


2-  I 3u(u'. 


(36) 


n M cos  (1)  t , - / , 

B = -2 ^ 3u  (u^ 

• 2 Ti  Y d® 


(37) 


B = 
X 


li  M cos  u)  t 
o o 

2 TT  Y 


[(2  - u^)  (u^  +1)"^'^^ 


(38) 


B = 
v 


Lt  M cos  U)  t 
o o 

o 2 ,4 

2 TT  Y d 


2 .,y-5/2 


(u^  +1) 


(39) 


wheiu 

X 

U = -^  (40) 

and  X and  d are  defined  in  Figure  III-l.  The  distance  variable  x , 

P P 

measured  from  the  location  of  the  moving  dipole  to  the  point  of  closest 
approach  to  the  sensors < is  assumed  to  vary  with  time  according  to  the 
expression 


X = -vt  (41) 

P 


A cosinusoidal  time  factor  has  been  included  in  Equations  (35)  through 

(39).  The  radian  angular  frequency,  cm  = 2 irf  , is  the  fundamental  of 

o o 

the  quasisinusoidal  emission  signals. 

For  the  range  computations  it  is  assumed  that  the  component 

signals  are  imbedded  in  background  noise.  Non-physically  realizable 

matched  filters  are  assumed  to  be  used  to  extract  the  signal  from  the  noise. 

A matched  filter  has  a frequency  response  that  is  the  complex  conjugate  of 

of  the  Fourier  transform  of  the  signal  to  which  it  is  matched.  In  order  to 

design  the  filters  necessary  to  extract  the  signals  of  Equations  (35)  through 

(39)  it  is  necessary  to  obtain  the  Fourier  Transforms  of  those  expressions. 

At  this  point,  it  is  noted  that  the  component  electric  field  signal  E may 

!! 

be  derived  from  B by  multiplying  by  y/p  • similar  fashion,  E may 
1 o X 

be  obtained  from  B by  multiplying  by  the  same  factor.  The  quantity 

II 

y/jj^  is  actually  l/p^  times  the  wave  impedance  at  the  surface  of  the  earth 
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looking  down.  It  can  bo  shown  that  for  noises  originating  in  sources 
sufficiently  far  away,  the  orthogonal  horizontal  components  of  the  electric 
and  magnetic  fields  are  related  in  exactly  the  same  way.  This  means  that 
the  ratio  of  the  signal  to  the  noise  is  the  same  for  the  electric  field  com- 
ponents as  it  is  for  those  of  the  magnetic  field.  It  is  redundant  to  derive 
the  filters  and  the  maximum  detection  ranges  for  the  components  of  both 
fiveldc.  For  this  reason,  the  calculations  are  continued  only  for  the  three 
components  of  the  magnetic  field . 

In  order  to  obtain  the  Fourier  Transforms  of  Equations  (37) 
through  (39)  it  is  noted  that  each  of  those  expressions  may  be  written  in 
the  form 

B(t)  = k h(t)  cos  u)  t (42) 

o 

where 

h(t)  = f[u(t)j  (43) 

and  f(u)  is  the  appropriate  function  of  u in  the  square  brackets . The 
factor  k in  (42)  includes  all  of  the  terms  outside  the  brackets  that  are 
constant  Vi^ith  respect  to  u and  t.  The  transform  of  (42)  may  be  obtained 
as 

B(uj)  = ~ [H(u)-«j  ) + H(ti)  + u)  )]  (44) 

Z o o 

where  the  transform  H(u>)  is  most  conveniently  written  in  terms  of  the 
Fourier  Transform  of  f(u)  as 
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In  (45) 


X = 


F(X)  = J f(u)  e ^ du 


is  the  expression  for  the  Fourier  Transform  of  f(u).  The  superscript  * 
with  F in  (45)  denotes  the  complex  conjugate. 

The  three  functions  of  u in  the  square  brackets  on  the  right 
hand  sides  of  Equations  (37)  through  (39)  are  designated  f^(u),  ^2^^^ 
and  f0(u),  respectively.  Thus 

fj(u)  = 3u(u^  + l)'^/^  {< 


f^(u)  = (2  - u^)  (u^  + 1)"^/^ 


fg(u)  = (U^  + 1) 


The  Fourier  Transforms  of  (48),  (49)  and  (50)  are  tabulated,  [Bateman,  1954], 


(51) 


Fj(X)  = -jX^Kj(X) 

F-(X)  = X^K  (X)  + XK,{X)  (52) 

Fg(X)  = 1/3  [ X^K^(X)  + 2XK^(X)|  (53) 

The  functions  F and  F and  the  magnitude  of  F.  are  shown  plotted 

M o X 

for  positive  X in  Figure  39.  It  should  be  noted  that  F2  and  F^  are  even 
functions  of  X v»hile  F^  is  odd. 

In  order  to  complete  the  Fourier  Transforms  of  the  three  com- 
ponents of  the  magnetic  field  it  is  necessary  to  first  substitute  (51)  through 
(53)  into  (45) . With  the  aid  of  (45)  the  results  are  substituted  into  (44) . 

Th“.  final  result  obtained  for  B is  written  as  an  example  of  the  set  of 

1 

three.  The  expression  is 


The  frequency  response  of  the  filter  matched  to  the  jl  component  of  the 
magnetic  field  is  the  complex  conjugate  of  B (uo).  Since  the  expression 

X 

of  (54)  is  real  the  equation  for  the  matched  filter  response  is  identical  to 
(54). 

The  response  of  the  filter  described  by  (54)  is  seen  to  be  the 
sum  of  two  filter  responses.  Each  has  a shape  defined  by  F^  but  one 
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filter  is  shifted  to  the  frequency  w while  the  other  is  shifted  to  - uj  . 

o o 

With  values  for  d,  v and  u)  typical  of  the  ship  detection  problem,  the 

o 

response  of  the  filter  translated  to  either  one  of  the  two  frequencies  has 
negligible  amplitude  at  the  other.  Thus,  the  shapes  of  the  matched  filters 
are  for  all  practical  purposes  the  same  as  those  described  by  the  functions 
Fj,  F2  and  F^.  The  term 

X = — (to  ± u)  ) (55) 

V o 

represents  a normalized  frequency  variable.  For  this  reason,  the  cur\^es 
of  Figure  39  can  be  interpreted  as  the  response  of  the  matched  filters  on 
one  side  of  the  center  frequency  and  the  functions  F^ , and  F^  are 
referred  to  as  the  matched  filter  shape  functions.  All  of  the  filters  may 
be  described  as  narrow  band  pass  filters  centered  on  the  fundamental 
frequency  of  the  quasisinusoidal  emission  signals 

The  output,  Sj,(t) , of  a filter  matched  to  a signal  S(t)  may 
be  written  as 


®F^^^  " f 1 e'^'^^dt  (56) 

— CO 

It  is  assumed  that  this  output  maximizes  at  t = 0 so  that  the  maximum 
output  of  the  filter  can  be  written  as 
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If  the  matched  filter  responds  to  noise  with  power  density  Pj^(u))/  the 
expected  value  of  the  rms  (standard  deviation)  of  the  output,  np(t) , is 
given  by 

03 

rms  = j Pjj((«)  I S(u))  1^  duD  (58) 

— 03 

The  detection  problem  amounts  to  comparing  the  character  of  noise  vvith 
rms  value  given  by  (58)  to  the  maximum  output  of  the  filter  produced  by 
the  signal  and  expressed  by  (57) . 

Two  simplifying  assumptions  are  used  for  the  derivations  to 
follow.  In  the  first,  the  power  density  Pj^(u))  is  assumed  to  be  a 
relatively  slowly  varying  function  of  frequency  similar  to  the  measured 
curves  shown  in  Figures  26  through  28.  The  signal  spectrum  S(a')» 
typical  of  the  component  emission  signals,  has  appreciable  amplitude  only 
over  two  very  narrow  bands  of  frequencies  centered  on  ± u)  . The  implica- 


tion is  that  Pjj(«))  is  essentially  constant  in  the  vicinity  of  the  frequencies 


where  Siw)  has  values  so  that  Equation  (58)  may  be  written  as  | 


The  remaining  calculations  arc  also  simplified  with  the  second  assumption 
that  the  filtered  noise  is  random  with  a Gaussian  amplitude  probability 
distribution. 

The  maximum  range  of  detection  problem  is  described  specifically 
as  follows.  The  probability  with  which  an  independent  point  of  np(t) 
is  not  to  exceed  the  value  of  Sp(0)  is  specified.  The  maximum  range  at 
which  this  requirement  is  satisfied  is  defined  to  be  the  maximum  range  of 
detection.  This  range  will,  of  course,  be  a function  of  the  specified 
probability . 

The  probability  specification  is  written  as 

Pr  |-^(0)<  np(t)  < +Sj.(0)|  = P (60) 

where  P is  the  specified  probability.  The  specification  is  chosen  to 
be  two  sided  for  the  following  reason.  The  filtered  signal  output  is 
quasisinusoidal,  even  more  so  than  the  unfiltered  signal.  Although  the 
maximum  of  the  filtered  signal  theoretically  occurs  as  a positive  peak  at 
t = 0,  in  practice  the  maximum  applies  to  a time  interval  representing 
several  cycles  of  the  sinusoidal  variations.  Thus,  a filter  output  that  is 
more  negative  than  -8^.(0)  is  considered  also  to  be  a significant  detection. 
Hence,  the  probability  specification  of  (60)  is  written  two  sided. 
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Dividing  each  of  the  terms  inside  the  braces  on  the  left  hand 


side  of  (60)  by  the  rms  value  of  the  noise  gives 


Pr 


np(t) 

-C  < — — < C 

rms 


= P 


where 


C 


Sp(0) 

rms 


(61) 


(62) 


The  independent  sample  points  of  the  functi-'n  (np(t)/rms)  have  zero 
mean  and  unit  variance  so  that  the  values  of  C for  a specified  value  of 
P can  be  obtained  from  tabulated  values  of  the  normal  probability  integral 
[Brunk,  1965].  Three  different  values  of  C are  considered  here.  They 
are  P = 0.9,  P = 0.99  and  P = 0.999.  The  corresponding  values  of 
C are  subscripted  with  the  complement  to  P expressed  as  a percentage. 


^10% 

1.645 

^1% 

= 

2.326 

(63) 

s.1% 

= 

3.291 

The  subscript  refers  to  the  probability  that  an  independent  sample  point  of 
np(t)  will  exceed  the  established  detection  threshold. 

Computation  of  the  maximum  detection  range  for  a given  value 
of  C begins  with  Equation  (62) . Substituting  for  Sp(0)  from  (57)  and 
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for  rms  from  (59)  and  then  rearranging  gives 


I S((jj)  1^ 


(64) 


The  subscript  p indicates  the  selected  value  of  C. 

Equation  (64)  is  now  evaluated  with  the  aid  of  Parseval's 
theorem  written  as 


00  00 

I S(u3)  d(i)  = J*  S^(t)  dt  (65) 

— 03 ' — CO 

The  right  hand  side  of  (65)  is  evaluated  for  S(t)  equal  to  each  of  the 

functions  B , B , and  B given  in  Equations  (37),  (38)  and  (39).  One 
li  1 V 

approximation  is  used  to  evaluate  the  resulting  integrals.  It  is  illustrated 
by  writing  the  integral  for  the  general  form  of  the  signal  function  given 
in  (42) . 


I = 


J 


— 00 


2 2 2 
K h (t)  cos 


U)  t dt 
o 


(66) 


The  envelope  function  h(t)  varies  slowly  with  time  compared  to  the 

oscillations  of  cos  w t.  In  particular  it  is  assumed  that  h(t)  can  be 

o 

considered  essentially  constant  over  any  one  half  cycle  of  the  cosine. 


4 

4 

ii 


I = J K^h^W  dt  « f cos^io^t  dt  (67) 


where 


t,  - t = — 

12  ()j 


(68) 


and  the  choice  of  t is  immaterial  as  long  as  it  is  between  t and  t . 

X Z 

The  integral  on  the  right  hand  side  of  (67)  can  be  evaluated  to  give 


2 2 

r h"'  (tg) 


(69) 


Since  h(t)  is  assumed  almost  constant  over  the  interval  (t  , t^) , 

X ^ 


h^(tg)  - tj)  « J h^  (t)  dt 

^1 


(70) 


so  that 


£ 

I 


h^  (t)  dt 


(71) 


If  the  whole  t axis  is  subdivided  into  consecutive  intervals  all  of  which 
have  equal  length  corresponding  to  one  half  cycle  of  the  cosine  then 


I « ^ 


J K^h^t)  dt 


(72) 
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The  substitutions  given  in  (40)  and  (41)  are  used  with  (72)  to  obtain  finally 


I S(ou)  du)  j f^(u)  du 


The  proper  forms  of  K and  f from  each  of  Equations  (37),  (38)  and  (39) 
are  used  in  (73)  and  the  resulting  integrals  evaluated  [Dwight,  I960]. 

The  results  are  substituted  into  (64)  to  obtain  the  three  detection  equations 


d = 
P 


/ M f / 45  X 10~^  Y 

\%)  V512TTaPj^ 


d = 
P 


/_M  ( 123  X 

W)  \512naPj^ 
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It  is  important  to  note  here  that  the  values  of  P^^  used  in  (74),  (75)  and 
C^6)  are  the  two  sided  power  density  functions  defined  for  both  positive  and 
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negative  frequencies.  Most  power  density  spectra  computed  from  ex- 
perimental data  including  those  shown  in  all  of  Figures  16  through  Figure  28 
are  one  sided  and  assume  that  all  of  the  power  is  in  the  positive  frequencies. 
It  is  therefore  necessary  to  divide  the  one  sided  power  density  spectrum 
values  by  Z if  those  values  are  to  be  used  in  Equations  (74),  (75)  or  (76). 

It  is  also  important  to  emphasize  the  units  used  for  the  power 
density.  MKS  units  are  required  for  Equations  (74),  (75)  and  (76).  Most 
experimentalists,  however,  follow  the  geophysical  convention  of  using  the 


unit  Y = 10  Wfebers  per  square  meter  for  the  magnetic  field  and  the  unit 
(mv/km)  for  the  electric  field.  These  geophysical  units  have  been  used 


in  Figures  16  through  28  so  that  the  appropriate  scaling  must  be  done 
before  any  values  are  entered  into  the  detection  equations. 

A review  of  the  three  Equations  (74),  (75)  and  (76)  shows  that 
B is  most  likely  to  give  the  largest  detection  range.  The  only  differences 

X 

between  (74)  and  (75)  are  the  terms  in  the  numerator  of  the  second  parenthesis 

of  each  equation  and  the  different  component  noise  power  densities.  In 

general,  there  is  no  reason  to  expect  that  the  noise  power  density  is 

different  in  either  axes.  This  means  that  the  larger  numerator  factor  of  (75) 

causes  that  equation  to  yield  a larger  value  for  d than  (74) . The  power 

P 

density  of  the  vertical  component  of  the  magnetic  field  noise  may,  however, 
be  quite  a bit  smaller  than  the  horizontal  components.  The  extent  to  which 


it  is  smaller  depends  upon  a number  of  factors  including  the  distribution  of 
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conductivity  in  the  solid  earth  beneath  the  sun  floor  and  the  depth  of  the 

overlying  water.  In  some  instances,  the  small  vertical  field  noises  could 

cause  (76)  to  give  a better  detection  ranca  than  (75).  However,  the 

vertical  field  signal  falls  off  more  rapidly  with  range  than  does  either  of 

the  horizontal  signals  so  it  appears  that  (75)  would  in  general  provide  the 

largest  value  of  d . 

P 

Assuming  this  to  be  the  case.  Equation  (75)  is  evaluated  with 

the  parameters  observed  for  the  merchant  ships  in  Port  Aransas  under  two 

sets  of  conditions.  The  first  is  the  set  that  existed  during  the  Poit  Aransas 

measurements.  A conductivity  of  one  mho  per  meter  is  used  for  o and 

-5  2 

P,,  has  the  value  10  v /Hz.  The  results  are  presented  in  Table  3. 

N 


The  second  set  of  conditions  is  considered  more  appropriate  to  the  open 

ocean.  The  velocity  of  each  ship  is  set  to  20  knots.  The  fundamental 

frequencies  of  the  emissions  are  scaled  from  the  value  measured  in  Port 

Amasas.  The  scaling  assumes  that  the  frequencies  are  proportional  to  the 

velocity  of  the  ship.  A conductivity  of  5 mhos  per  meter  is  used  to  represent 

-5  2 

that  of  the  sea.  water.  As  before,  a noise  power  density  of  10  y /Hz 

is  used  forP-.  . The  resulting  d values  for  the  different  merchant  ships 
N p 

i 

are  tabulated  in  Table  4. 

The  extent  to  which  the  noises  caused  by  the  surface  water 
waves  affect  the  detection  ranges  can  be  determined  from  the  curves  of 
Figures  32  through  38.  Consider  specifically  the  case  of  a magnetic  field 
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measurement  made  near  the  surface.  The  curves  of  Figure  33  show  Die  power 
density  spectra  of  magnetic  field  noises  computed  for  the  spectral  distribution 
of  wave  height  energy  shown  in  the  curves  of  Figure  32.  The  parameter  in 
both  plots  is  wind  speed.  A+  frequencies  above  IHz  the  wave  induced 

magnetic  field  noise  power  density  is  essentially  independent  of  wind  speed. 

-5  2 

The  level  is  below  the  value  of  10  y /Hz  used  for  the  detection  range 
calculations  leading  to  the  results  of  Tables  3 and  4.  This  means  that  the 
wave  induced  noises  should  not  materially  reduce  the  maximum  detection 
ranges. 

-5  2 

It  should  be  pointed  out,  however,  that  a value  of  10  y /Hz 
is  a typical  summertime  low  latitude  value  and  does  not  represent  a minimum. 
In  some  instances,  for  example,  during  the  winter  months  in  polar  regions, 
the  power  density  of  the  background  noises  may  be  considerably  less.  In 
this  event,  the  ocean  wave  induced  noises  may  indeed  serve  to  limit  the 
maximum  detection  ranges.  The  detection  range  can,  however,  be  increased 
by  measuring  tJ  ,i  agnetic  field  at  some  depth  below  the  water's  surface. 

The  effect  can  be  noted  with  the  aid  of  Figure  36.  That  figure  shows  curves 
of  the  ocean  wave  induced  magnetic  field  power  density  spectra  at  a depth 
of  100  meters.  The  level  of  power  density  at  frequencies  in  the  1 to  10  Hz 
range  is  seen  to  be  negligible  at  this  depth. 

The  signals  and  background  noises  also  attenuate  with  depth 
in  the  sea  water  but  nowhere  near  as  rapidly  as  the  wave  induced  noises. 
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A measure  of  the  signal  attenuation  is  obtained  by  noting  that  a skin  depth 
at  5 Hz  in  sea  water  is  roughly  100  meters.  This  means  that  at  5 Hz 
the  signal  power  at  100  meters  would  be  reduced  by  a factor  of  about  7. 
Frequencies  below  5 Hz  would  be  attenuated  to  a lesser  extent.  It  is 
concluded  from  the  above  that  for  a given  set  of  wave  conditions  there  is 
some  depth  of  water  below  which  the  wave  induced  noises  would  have 
negligible  effect  on  the  detection  ranges.  The  limiting  factor  would  be 
the  level  of  natural  background  noises. 
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VI.  CONCLUSIONS 

This  study  has  been  concerned  with  both  the  theoretical  and 
experimental  aspects  of  detecting  the  electromagnetic  emissions  from  ships 
in  motion.  Measurements  made  on  merchant  ships  in  the  Corpus  Christi  ship 
channel  clearly  indicate  that  the  emission  signals  behave  like  the  fields 
of  a horizontal  electric  dipole  at  the  position  of  the  ship. 

Based  upon  this  fact,  detection  ranges  are  computed  for  the 
case  where  the  signals  are  extracted  from  the  background  noises  with  a 
matched  filter.  The  ranges  obtained  depend  upon  the  effective  strength 
of  the  moment  of  the  source  of  the  omissions.  This,  of  course,  varies 
from  ship  to  ship  and  may  vary  with  time  on  a given  ship  as  the  condition 
of  the  hull  and  the  screws  changes.  Several  of  the  ships  passing  through 

the  ship  channel  produced  no  mea';urable  emission  signals.  The  assumption 
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that  they  could  not  be  detected  may  not  bo  exactly  correct.  If  the 
measurements  had  been  processed  with  a sol  of  matched  f’  tors  there  may 
have  been  a low  level  of  emission  that  could  have  been  detected.  However, 
the  maximum  detection  ranges  computed  from  these  ships  under  typical  sea 
conditions  would  in  any  event  be  less  than  the  smallest  range  shown  in 
Table  4.  Such  a range  would  be  too  small  to  be  of  interest.  This  would 
seem  to  cast  some  doubt  on  the  efficacy  of  detection  methods  based  solely 
on  the  propeller  modulated  emissions.  They  should  not  be  completely 
discounted,  however,  because  the  detection  ranges  obtained  for  ships  with 
large  source  moments  have  quite  reasonable  values  at  around  2 kilometers. 

In  addition  it  is  noted  hero  that  all  of  the  detection  range 
calculations  are  based  on  single  point  measurements  of  the  component  fields 
associated  with  the  emissions.  In  previous  studios  [Fowler,  1973],  our 
laboratory  has  used  linear  filtering  techniques  with  measurements  of  the 
spatial  field  gradients  to  significantly  reduce  the  background  noise  levels. 
This  work  was  done  in  connection  with  studies  of  the  maximum  detection 
ranges  that  could  be  achieved  by  sensing  the  static  dipole  moments  of 
moving  ships.  The  frequencies  involved  were  well  below  those  of  the 
present  study.  Reductions  of  the  noise  levels  by  factors  of  up  to  40  db 
were  produced  with  filters  that  corrected  for  the  inhomogeneous  distribution 
of  conductivity  in  the  earth.  These  filters,  termed  "geology  filters", 
were  originally  developed  at  our  laboratory. 


Similiar  techniques  could  be  used  to  reduce  the  background  noise 
levels  at  the  frequencies  of  the  propeller  emission  signals.  A reduction 
of  the  noise  power  by  a factor  of  40  db  would  increase  the  detection  ranges 
listed  in  Tables  3 and  4 by  the  one  fifth  power  of  that  factor.  This  would 
amount  to  an  increase  of  6.3  times  and  would  produce  a very  usable  detection 
range  under  sea  conditions  of  over  10  km  for  the  ships  with  large  source 
moments.  The  smallest  detection  range  listed  in  Table  4 would  increase 
to  a very  respectable  value  of  over  3 krn. 

With  this  projected  result  in  mind,  the  aithors  recommend  that 
the  feasibility  of  implementing  such  a system  be  investigated.  On  the 
basis  of  results  of  previous  studies  it  is  also  recommended  that  such  a 
system  incorporate  the  additional  features  necessary  to  detect  the  static 
dipole  moment.  This  could  be  done  with  the  same  set  of  sensors. 
Simultaneous  processing  to  extract  both  the  static  and  propeller  induced 
moments  should  produce  a detection  system  with  significant  capability. 
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TABT.K  1 

SUMMARY  or  DATA  RUNS 


Date 

Description  of  Run 

Run  No. 

09-18 

Giovanella  d'Amico 

091873-1 

Wanderer  & David  P.  Reynolds 

091873-2 

Honner 

091873-3 

Calibrate 

091873-4 

09-19 

Ero 

091973-1 

Melias 

091973-2 

Calibrate 

091973-3 

Stainless  Trader 

091973-4 

David  P.  Reynolds 

091973-5 

09-21 

Benja  River 

092173-1 

Pitria  Star 

092173-2 

09-24 

Background 

092473-1 

Background  (No  Digital  Tape) 

092473-2 

Chemical  Venturer 

092473-3 

09-25 

Background 

092573-1 

Background 

092573-2 

Cali  brate 

092573-3 

Inger 

092573-4 

Ero 

092573-5 

Tumi 

092573-6 

Richard 

092573-7 

Background 

092573-8 

09-26 

Texas  Sun 

092673-1 

Background 

092673-2 

Background 

092673-3 

Richard 

092673-4 

n 3 r'V/nrrMinH 

092673-5 

David  P.  Reynolds 

092673-6 

Background 

092673-7 

Background 

092673-8 

Aghios  Spyridon 

092673-9 
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TABLE  1 
(continued) 


Date 

Description  of  Run 

Run  No. 

09-28 

Mariotte 

092873-1 

Background 

092873-2 

Background  (Mo  Digital  Tape) 

092873-3 

Background 

092873-4 

Hellenic  Ideal 

092873-5 

09-30 

Master  Stefanos 

093073-1 

Marjorie  Lykes 

093073-2 

Solholt 

093073-3 

10-02 

La  Bonita 

100273-1 

Wingull 

100273-2 

Gulfcr*  St 

100273-3 

Master  Stefanos  & Aquacharm 

100273-4 

10-03 

Calibrate 

100373-1 

Nopal  Tellus 

100373-2 

Monmouth 

100373-3 

10-04 

Electrode  Test 

100473-1 

Hess  Trader 

100473-2 

David  P.  Reynolds  & Coastal  Texas 

100473-3 

10-05 

Richard 

100573-1 

Texaco  Oregon 

100573-2 

10-06 

Ship  (Name  Not  Recorded) 

100673-1 

10-07 

Rio  Balsas 

100773-1 

Atlantic  Prestige 

100773-2 

Nopal  Trader 

100773-3 
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TABLE  3 


TABLE  4 

DETECTION  RANGES  - MATCHED  FILTER 
SEA  CONDITIONS 
id  = 5 v/m,  V = 20  Knots, 


d 10% 


d 1% 


d 0.1% 
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fon 


Figure  8.  Recorded  Time  Traces  - Texas  Sun 


Figure  10.  Recorded  Time  Traces  - Mariotte 
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POWER  SPECTRA  - BACKGROUND  RUN 
Figure  26.  Background  Noise  092673-002 

71 


o i 


* ^ •;  ^ 

w o 

II  ‘ 

_ E 

^o 


UJ 


C)  . 


I ^-1 
0.  „ ^ 


•Ji 

o . 


10 


-3 


RUN  NO  092873-002 


o « 


O O 


O O 


T r-T*i  -i-rf-rr 


10^? 

POWER  SPEClRft 


io-> 

FRLQUINCT 


TT-I-I ! l-T  -t-1 

10° 


BRCKGROUK'D  RjN 


-'3 

C>i  4 

R-i 

o a 

2l 


•8 

o » o 


RUN  NO  092873-002 


£1 


o.'' 

XoT  1 

£ I 

S o.' 

2 

O. 

iir 


o.|. 

10-  3 


O ■ « 

* <5  o I 


1 rj  r »— »— 7-T  T rr?— — t — v 

!0'^  ]0-' 

f R(  autNC7 


r ‘**T — I — r-T  • - 1 1 


POWER  S^EETRO 


10° 

Rm^KGROUNI;  Rl!\ 


Figure  28.  Background  Noise  092873-002 
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Appendix 

Wave  Noise  Derivation 

Expressions  for  the  electric  and  magnetic  fields  generated  by  ocean 
wave  motions  of  the  conductive  sea  water  perpendicular  to  the  geomagnetic 
field  are  derived  in  this  appendix.  The  derivation  assumes  that  the  ocean 
wave  motions  are  described  by  the  simple  Airy  theory  of  surface  waves. 

No  attempt  has  been  made  to  include  the  effects  of  internal  waves.  For 
other  treatments  of  the  electromagnetic  fields  generated  by  ocean  waves, 
the  reader  is  referred  to  references  f Weaver,  1965],  fBeal  and  Weaver, 
1970],  [Podney,  1974A],  fPodney,  197  iB]. 

Simple  Airy  wave  theory  assumes  that  the  air  sea  interface  varies 
sinusoidally  about  mean  sea  level  with  respect  to  both  time  and  horizontal 
distance.  If  the  position  of  the  surface  with  respect  to  mean  sea  level  is 
designated  by  r\ , then 

/ .,1  " kx) 

Ti(x,t)  = ae"' 

This  represents  a wave  propagating  in  the  +x  direction  whose  motions  are 
independent  of  y as  illustrated  in  Figure  lA,  If  v represents  the  velocity 
of  the  water  at  any  point  below  the  surface,  we  may  write 

v^a  u+dw 

The  Airy  wave  theory  ossumes  that  the  flow  of  water  associated 
with  the  wave  is  irrotational  expressed  by 


Vx  V = 0 
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and  that  the  water  is  incompressible,  expressed  by 


V.  v=  0 A-1 

The  assumption  of  irrotational  flow  suggests  the  use  of  a scalar  velocity 
potential,  , given  by 


v=  -V'|f  A-2 

Substituting  A-2  into  A-1  gives  La  Place’s  equation, 

^2  1.  = 0 A-3 


for  the  velocity  potential. 

Solutions  of  A-3  may  be  written  as 

L /A  kz  , „ -kz.  j(ujt  - kx)  - . 

ilf  = (Ae  + Be  ) o . A-4 

In  order  to  evaluate  A and  B we  use  the  kno'Am  values  of  the  vertical 

components  of  the  water  velocity,  w,  at  z = 0 and  at  z = -h.  In  terms 

of  , w may  be  expressed  as 

1 /A  „ -kz.  j((t't-  kx) 
w=-~ — = -k(Ae  -Be  )e 
az 


At  z = 0 , 


ar|  . j{t«t  - kx) 
w = = juua  e 

ot 


so  that 


(A  - B)  = -Jf 


AL  z = -h,  0 giving 


. -kh  n kh 
Ae  -Be  =0 


A- 5 


A-6 
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Solving  A-  5 and  A-6  for  A and  B gives 


A = - 


2k  sinh  kh 


P ^ e 

2k  s’.ih  kh 

Using  A-7  and  A-8  in  A-4  we  obtain 

CO  sh(2H-  h)  j(a>t-  kx) 

’ k sinh  kh  ^ 

Substituting  this  expression  for  in  A-2  yields  for  the  component 
velocities 


cosh  k(z+h)  )(u)t  - kx) 
u = u)a — T-rTT — ^ e 
sinh  kh 


sinh  k(z+h)  ){u!t  - kx) 
v=juia — r-rTT — e 
Sinh  kh 


In  order  to  obtain  an  expression  for  k in  terms  of  x (the  dispersion 
relation)  we  make  u"e  of  Euler's  equation  (equation  of  motion).  Thus 
dv 

P^=  - pvu  - Vp 


where  pis  the  mass  density  of  the  water. 


u = gr  t-  C, 


is  the  gravitational  potential  with  g 9.81  m/sec  and  Cj  a constant; 
and  p is  the  scalar  water  pressure.  Expanding  the  total  derivative  gives 


p(-^  + V*  Vv)  = - pVu-Vp 
or 
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We  linearize  this  equation  by  assuming  that  the  waves  have  a reasonably 
small  amplitude  so  that  the  factor  v • Vv  is  negligible.  There  remains  the 


equation 


d’Jf 


- pv  — = - oVu-Vp 
ot 


that  may  be  integrated  to  obtain 
S ’i( 


M --  pgz  + c^-p+c^ 


A-10 


where  is  a constant  with  respect  to  the  space  coordinates. 


At  the  water's  surface  z - r\  and  p = p is  the  sea  level 

a 


atmospheric  pressure,  assumed  to  be  constant  with  respect  the  space 
coordinates.  Evaluating  A-10  at  the  surface  gives 


9t 


lz=„=  - >1+  <=1  - Pa  + S 


Since  and  p are  functions  of  x and  C, , p and  C„  are  not,  it  must  be 

la  z 


true  that 


Cj  - + C3  = 0 


and 


Ji.  I 

at  'z=Ti 


A-11 


We  can  again  use  the  assumption  of  reasonably  small  amplitude  waves 
to  advantage  by  noting  that  for  such  waves, 

^ (x,  n(x)  ) (x,o)  A-12 
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Substituting  A-12  into  A-11  gives 


11 

ht 


gn 


Substituting  A-9  into  A-12  gives  for  the  dispersion  relation 


k 


coth  kh  = g. 


A-13 


At  this  point  we  use  an  assumption  that  results  in  considerable 
simplification  of  the  development  to  follow.  We  assume  that  the  water 
is  of  sufficient  depth  (h  large  enough)  so  that 


coth  kh  1 


and 


A-14 


Also  letting  h be  very  large  in  A-9  obtains  for  the  velocity  potential 


jj'a  kz  j(u)t  - kx) 

, _ g g 

k 


A-15 


Equations  A-14  and  A-15  describe  the  behavior  of  waves  where  the  water 
is  deep  enough  to  attenuate  the  effect  of  waves  reflected  from  the  ocean 
floor. 


The  interaction  of  the  moving  sea  wab  . ■ -ing  conductivity  a with 

the  constant  geomagnetic  field  B^is  described  oy  the  magnetohydronamic 
form  of  Maxwell's  relations 

V X E = - juup oH  A-16 


4 

;| 

*5 

1 

I 

i3 

a 


$ 

3 

I 

$ 

■j 

u 

ft 

n 

s 

I 
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A-17 


Vx  H = J 

f = 0 (E  + V X B)  A- 18 

where  H,  E,  and  J are  the  magnetic  field,  the  electric  field  and  current 
density  respectively.  The  quantity  is  the  magnetic  permeability  of 
free  space.  In  equation  A-17  we  have  assumed  that  the  displacement 
currents  are  negligible  and  this  is  quite  justified  at  the  frequencies  and 
conductivities  involved  in  this  study. 

Taking  the  divergence  of  equation  A-17  we  obtain 

V-  J = V . (Vx  H ) - 0 A-19 

Taking  the  divergence  of  A-18  and  using  the  result  in  A-19  gives 

V.  J=ov.  (E+vxBj  s 0 A-20 

Here  we  have  assumed  that  o is  a constant  and  in  this  application  this 
appears  reasonably  valid.  Using  A-2  for  v in  the  term  v x B^  in  A-20 
gives 

vx6  =-Vi!(xB  =-VxilfB  A-21 

o o o 

where  the  right  hand  equality  results  from  a vector  identity.  Solving 
for  V . E from  A-20  and  substituting  from  A-21  yields 

V.E=-V*  (vxB  ) = -V.  (-VxiliB  ) =0 
o o 

E is  therefore  solenoidal  and  may  be  expressed  as  the  curl  of  a vector 
- £ 

potential,  A . Thus 
E=  Vx  A^ 
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A-22 


Taking  the  curl  of  A-16,  substituting  for  Vx  H from  A-17  and  for  E from 
A-22  and  substituting  for  v x from  A-21  gives 


_E  -E 

Vx  Vx  Vx  A +jO)pavxA  = juDpa  (Vx  * B ) 

o 


A- 23 


Integrating  A-23  and  using  a vector  identity  yields 

E -E  2 -E  - 

7$+v(^*A  )-V  A + ju)(icrA  = juu|a'!iB^o 

where  $ is  an  arbitrary  function  of  the  space  coordinates  that  may  be 
chosen  for  convenience.  To  this  end  we  choose 
$ = - V . A 


to  obtain 


7^  A^  - ju)iaaA^= B a A-24 

o 

Equation  A-24  is  an  inhomogeneous  equation  for  the  ’ *or  electric 

_E 

potential  A with  a source  term  that  involves  the  motion  of  the  conductive 
sea  water  interacting  with  the  geomagnetic  field. 

Rather  than  attempt  to  derive  all  of  the  electromagnetic  fields  from 
the  vector  electric  potential  we  will  use  that  potential  to  derive  only  a 
portion  of  them.  In  particular  we  will  obtain  from  the  vector  electric 
potential  only  those  fields  that  can  be  derived  from  a solution  of  A-24 
that  has  only  a y component.  Reference  to  A-24  shows  that  this  means 
that  we  are  seeking  only  those  fields  excited  by  the  y component  of  the 
geomagnetic  lield.  The  fields  excited  bv  the  x and  z components  of  the 
. omagnetic  field  are  more  conveniently  derived  from  the  vector  magnetic 
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potential  as  will  be  discussed  later.  For  the  present  case  we  assume 
A^(x.z)  = Ay(x,z) 

_ £ 

is  a solution  to  A-24,  The  electric  field  may  be  obtained  from  A with 

the  aid  of  equation  A-22.  An  expression  for  the  magnetic  field  in  terms 
- E 

of  A may  be  obtained  by  first  solving  equation  A-16  for  H and  then 
substituting  for  E from  A-22.  This  yields 

H = - V X vx  A^  = - (v  (v  • A^)  - A^)  A-25 

JUUli  J{10^ 

where  the  right  hand  equality  results  from  an  identity  for  the  vector 

differential  operators.  Further  simplification  results  if  we  make  use  of 
~E 

the  fact  that  A has  only  a y component  that  is  independent  of  y.  This 

implies  that 
-E 

V • A = = 0 

By 

so  that 

E=  V^A^  A-26 

JUUji 

2 -E 

Substituting  for  V A from  A-24  gives 

H = a (A^  - Ilf  B a ) A-27 

oy  y' 

This  simple  expression  shows  that  the  magnetic  field,  H,  derived  in  this 

manner  has  only  a y cc>mponent.  Since  the  y component  is  transverse  to 

the  direction  of  propagation  of  the  ocean  waves  and  to  the  water  motions 

- - -E 

v.ltnin  tne  wave  we  will  refer  to  the  E and  H fields  derived  from  A as 


Transverse  Magnetic  (TM)  fields. 


In  order  to  derive  the  rest  of  the  electromagnetic  fields  excited  by 


the  ocean  wave  motions,  we  note  that  the  magnetic  field  is  solenoidal, 
as  may  be  verified  by  taking  the  divergence  of  A-16,  to  write 

B = H H = DxA  A-28 

o 

where  A is  vector  magnetic  potential.  Substituting  the  expression  for  H 

obtained  from  A-28  into  A- 17  gives 

7 X Vx  A = M J A-29 

o 

Substituting  for  J from  A-18  and  using  an  identity  for  the  vector  differential 

operations  results  in  the  equation 

V (7  ' A ) - A = li  a (E  + V X B) . A-30 

o 

We  obtain  an  expression  for  E in  terms  of  A by  substituting  A-28  into 
A-16  to  get 

7x  E = - juu7x  A 

Integrating  both  sides  results  in 

E = - ju)A  - 7$  A-31 


where  f is  an  arbitrary  function  of  the  space  coordinates  that  may  be 
chosen  for  convenience.  We  choose 

7 • A = - no  $ A-32 


Substituting  this  and  the  expression  for  E from  A-31  into  A-30  we 
obtain 

2 - - - - 
7 A - jiDpioA  = - pia  V X B 


A-33 
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Here  as  before  we  assume  that  the  vector  magnetic  potential  has  a y 
component  only . Thus 

A (x,z)  = a^A^(x,2), 

where  A^{x,z)  is  a solution  to  the  y component  of  equation  A-33. 

Because  A is  independent  of  y 
^ 9A 

7 • A(x,z)  = ■ ^ = 0 

According  to  A- 32 
$=  0 


so  that  A-31  becomes 

E = - j(uA  A-34 

This  equation  shows  that  the  electric  field  derived  in  this  manner  has  only 

a y component.  Following  the  same  reasoning  used  to  identify  the  fields 

derived  from  the  vector  electric  potential  as  TM,  the  fields  derived  from 

the  vector  magnetic  potential  are  termed  Transverse  Electric  (TE) . 

We  summarize  the  derivations  to  this  point  with  a review  of  the 

equations  describing  the  electromagnetic  fields  generated  by  ocean  wave 

motions.  In  this  development  the  total  electromagnetic  field  is  separated 

into  TM  and  a TE  parts  that  are  derived  from  the  y components  of  the 
-E  - 

vector  potentials  A and  A respectively.  The  y components  of  the 
potentials  are  solutions  to  the  y components  of  equations  A-24  and  A-33, 
rewritten  here  in  component  form  as 


2 E E 

7 Ay  - ju)|aoAy  = - jouijc  ? 


A-35 
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and 


V^A  - juJiacrA  = -na(—  B - — — B ). 
y y Bx  oz  S z ox 


A-36 


The  velocity  potential,  , subject  to  the  deep  water  wave 


assumption  is  given  by  equation  A-15  restated  he;e  as 

, juja  kz  j('X’l  - kx) 

--i- e e 

k 


A-15 


The  TM  electric  and  magnetic  fields  are  obtained  from  the  vector 
electric  potential  with  equations  A-22  and  A-27  respectively.  The  TE 
electric  and  magnetic  fields  are  obtained  from  equations  A-28  and  A-34 
respectively.  We  now  proceed  to  obtain  solutions  to  equations  A-35 
and  A-36. 

Consider  first,  solutions  to  equation  A-35.  We  assume  solutions 
that  vary  with  x and  t in  the  same  manner  as  ’J*  as  indicated  in  A-15. 
Thus 


A-37 


Substituting  A-37  and  A-15  into  A-35  and  performing  the  derivative  with 


respect  to  x gives 

2 E 

^ - (k  + juuna)  A (z)  = - jx'pa’Jr  (z)  B 


dz^ 


oy 


A-38 


where 

(z)  = - e^^  A- 39 

2 

The  partial  second  derivative  with  respect  to  z implied  by  the  V operator 
in  A-35  has  been  replaced  by  the  total  second  derivative  in  A-38  since 
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Ay(z)  varies  only  'vith  z.  Equation  A-38  may  be  written  as 
dV(z)  + (2)  = f(z) 


A- 40 


where 


^2  , 2 . 

C = - k - jajfiCf 


A-41 


f(z)  = - jiDiJ  0 B 


A-42 


A particular  solution  to  A-40  may  be  written  as 


A (z)  = fG(zj  z ) f(z  ) dz 
py  '000 


A-43 


where  the  integral  is  performed  over  all  z where  f(z  ) is  non  zero  and 

o o 


where 


_i_  -iC  (z  - z ) 
r ~ ior  ^ o»2>ii 

G(z|z  ) = ■[  ...  . 


A-44 


Substituting  into  A-43  from  A-42  and  A-44  and  rearranging  terms  gives 


A®  (z)=-?L!iHA  B r^-iC2  r%(k+JC)2od, 
py  j2Ck  oyL®  o 


A-45 


J O J 


The  integrals  in  A-45  are  easily  evaluated  and  after  some  algebraic 
manipulation  we  obtain  for  the  particular  solution 


py  k 


Bpy  ^k+i£  ^jC  2 _ gkZj 


A-46 
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To  thio  solution  we  add  a complementary  solution  of  the  homogeneous 
form  of  equation  A- 40, 


^ (z)  = 0. 


A- 47 


in  order  to  match  the  electiomagnetic,  boundary  conditions.  Solutions  of 


A-47  may  be  written  as 
Asy  (z)  = G e^^  ^ ^ 


A-48 


where  G is  an  undetermined  constant  whose  value  is  to  be  chosen  to 
match  the  boundary  conditions.  Consider  now  the  bou  idary  conditions 
for  the  TM  fields. 

In  the  air  above  the  water's  surface  (z-'  O)  the  electrical  conductivity 
a is  assumed  to  be  zero  so  that  there  are  no  electrical  currents  and 
equation  A- 17  written  for  the  magnetic  field  intensity  in  the  air  (designated 


by  the  subscript  a)  becomes 


7x  H = 0 
a 


Equation  A-26  is  valid  in  the  air  space  for  H , 
a j(up  a 


But,  according,  to  A-24  if  a = 0 

=0 

a 

so  we  conclude  that 


H = 0 . 
a 


A-49 
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This  requires  that  all  components  of  H including  the  components  tangent 

a 

to  the  air  sea  interface  on  the  air  side  of  the  boundary  be  zero.  Because 

the  tangential  components  of  the  magnetic  field  intensity  are  continuous 

across  a boundary  we  conclude  that  on  the  water  side  of  the  interface  the 

components  of  H (as  considered  hero  there  is  only  a y component  of  the 

TM  fields)  tangent  to  the  interface  are  also  zero.  Using  equation  A-27  to 

express  the  magnetic  field  we  oL.ain  the  magnetic  field  boundary  condition 

in  terms  of  the  vector  electric  potential  as 

A ^ (o)  - vi' (o)  B =0  A- 50 

y oy 

£ 

where  the  x and  t dependence  of  A^  and  '!>  have  been  suppressed. 

We  now  add  the  particular  and  complementary  solutions  given  by 
equations  A-46  and  A-48,  substitute  the  result  into  A-50  and  solve  for 
G to  obtain 


oua  (k  +jC) 

^oy 


A-51 


Using  this  value  for  G we  write  the  expression  for  the  vector  electric 


potential  from  wliich  we  derive  the  total  electric  vector  potential  as 

r. 


Ay  (x,z,t)  = -■ 


_ kz  j(aot  - kx) 

Bee 

k oy 


A- 52 


Using  equations  A -22  and  A-27  to  obtain  the  electric  and  magnetic  fields 


gives 


E = jeua  B e^^' 

X — oy 

T,  T. 

E = - 'jua  B e 
z oy 
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A-53 

A-54 


H = 0 A-55 

y 

TE  Fields 

In  the  same  fashion  as  used  w?ith  the  TM  fields  we  assume  solutions 
of  A-36  that  vary  with  x and  t in  the  same  manner  as  \ii  as  indicated  in  A-15. 
Thus 

A (x,z,t)  = A (z)  A-5G 

y y 


Substituting  A- 56  into  A-36  and  performing  the  derivative  with  respect  to 
X gives 


d\(z)  2 

+ C Ay(z)  = g(z) 


dz^ 


As  before  C is  given  by  A-41  and 


A- 57 


g(z)  = pak  ’I'  (B^  + jB^ J 
OX  oz 


A-58 


A particular  solution  to  A- 5 7 may  be  written  as 

A (z)  = I* G(z I z ) g (z  ) dz 
py  V ' ' o'  ^ ' o'  o 


A- 59 


where  the  integral  is  performed  over  all  z^  for  which  g (z^  is  non  zero. 
The  quantity  G(zl  z^)  is  given  by  A-44. 

Substituting  into  A- 59  from  A-44  and  A-58  and  rearranging  terms 


gives 


A_,.(Z)  ..  (B._  - IB., ) [ ^ I 0'"+  dz 


py 


2C 


oz  ox 


A-60 


o dz  'I 
o ■' 
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The  integrals  in  A-60  are  easily  evaluated  and  after  some  manipulation 
we  obtain  for  the  particular  solution 

A'  + i r i \r  ’7 

A-61 


A (z)  = a(P  +JP  ) "] 

ny  — ox  oz  ]2C 


To  this  solution  we  add  a complementary  solution  of  the  homogeneous  form 
of  equation  A-57 

d\(z)  2 

+ Ckiz)  = 0 

dz^  y 


in  order  to  match  the  electromagnetic  boundary  conditions.  Solutions  of 

the  above  equation  may  be  written  as 

A (z)  = F e^^  ^ 
sy 

where  F is  an  undetermined  constant  hose  value  is  to  be  chosen  to  match 
the  boundary  conditions.  Consider  now  the  boundary  conditions  for  the 
TE  fields. 

In  the  air  space  above  the  water  the  y component  of  the  vector 


magnetic  potential  satisfies  equation  A-36  with  0=0.  Thus 


V A = 0 
ay 


A-62 


We  assume  solutions  to  A-62  that  vary  with  x and  t in  the  same  manner  as 
we  have  assumed  for  the  potentials  in  the  water  as  expressed  by  A-56. 
Substituting  into  A-62  gives 


- >=  - 0 


A-63 
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A solution  to  A-63  that  behaves  according  to 


Lim  A (z)  = 0 
ay 

Z 00 

is  expressed  by 

A (z)  = T e"^'^  A- 64 

ay 

where  T is  an  undetermined  constant. 

Both  above  and  below  the  water's  surface  the  horizontal  components 

of  the  magnetic  and  electric  fields  may  be  derived  from  the  appropriate 

vector  magnetic  potential  with  the  aid  of  equations  A-28  and  A-34 

respectively.  In  both  the  air  and  the  water  the  vector  potential  has  only 

a y component  so  that  equations  A-28  and  A-34  become 

H(a)x  - — ■ A-65 

Po 

and 

E(a)y  - juuA(a)y  (z)  A-66 

where  the  x and  t dependence  has  been  supressed.  The  subscript  (a)  is 

intended  to  mean  that  equations  A-65  and  A-66  apply  to  A^  in  the  water 

as  well  as  A ■ in  the  air.  In  order  to  evaluate  the  boundary  conditions 
ay 

for  the  horizontal  components  of  the  fields  just  below  the  water's  surface 
we  first  evaluate  those  fields  in  the  air  space  just  above  the  surface. 
Substituting  A-64  into  A-65  and  A-66  and  evaluating  the  result  at  z = 0 
gives 

H = — kT  A-67 

ax  p 

o 
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and 


E = - j(JuT. 
ay 


A-68 


Because  of  the  continuity  of  the  tangential  magnetic  and  electric  fields 


across  an  interface  we  may  also  write 


1 


kT 


A- 69 


and 


E = - ju)  T 

y 


A-70 


for  the  fields  on  the  water  side  of  the  interface.  Substituting  for  T in 
A-69  its  value  derr  ed  from  A-70  gives 


H = - : 

X 


- E 

y 


Substituting  for  and  from  A-65  and  A-66  and  rearranging  gives 


k A (o)  + I =0 

y 8z  'z=0 


A-71 


as  the  electromagnetic  boundary  condition  expressed  for  the  vector 
magnetic  potential.  If  we  now  add  the  particular  and  complementary 
solutions  A-51  and  A-63  and  substitute  the  result  into  A-71  with  z = 0 
we  obtain  an  equation  that  may  be  solved  for  G to  give 


G = 


(B  +jB  ) 


(k  - KV 


2jC  ox  ox  k + jC 


Using  this  value  for  G the  total  vector  magnetic  potential  written  as  the 


sum  of  the  particular  and  complementary  solutions  is  given  by 

hi  /n  . n \ r kz  2k  jCz,  j(t«t-kx) 

A (x,z,t)  = a (B  + jB  ) [e  — ro.  e ] e 

y ' - ox  ■’  oz'  *•  (k+jC) 


A-72 
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Using  equations  A-28  and  A-34  to  obtain  the  electric  and  magnetic  fields 
gives 


E = - j'xa  (B  + jB  ) [e^^ 
y - ox  oz  k+  jC 


H = - k(B  + jB  ) [e^^  - 

X H ox  ' oz'  '■  (k+  jC) 


. - . . r kz  2k  C z, 

H =-j k(B  + jB  ) [e  -,r- — : . e J 

z M ox  oz  (k+  jC ) 


A-73 

A-74 

A-75 


Equations  A-73,  A-74  and  A-75  give  the  TE  fields  generated  by  the  ocean 
wave  motions  perpendicular  to  the  geomagnetic  field  while  equations  A-53 


A-54  and  A-55,  previously  derived  give  the  TM  fields. 
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A series  of  measurements  of  the  EM  emissions  from  moving  merchant 
ships  were  made.  Analysis  of  the  recorded  signals  shows  that  they  correspond 
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the  ship.  Also  recorded  were  the  levels  of  the  natural  EM  background.  Using 
these  levels  the  maximum  ranges  at  which  a matched  filter  could  extract  the 
ship  emission  signals  are  derived.  Also  considered  are  the  effect  of  electro- 
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